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Abstract 

We calculate the optical and Raman response within a phenomenological model of fermion quasiparticles coupled to 
nearly critical collective modes. We find that, whereas critical scaling properties might be masked in optical spectra 
due to charge conservation, distinct critical signatures of charge and spin fluctuations can be detected in Raman 
spectra exploiting specific symmetry properties. We compare our results with recent experiments on the cuprates. 
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1. Introduction 

It has been proposed that the anomalous physi- 
cal properties of the cuprates arise from the proxim- 
ity to a quantum critical point (QCP) located near 
optimal doping [1,2,3,4]. The onset of the ordered 
phase, however, can be prevented by the low dimen- 
sionality of the Cu02 layers, by disorder, and/or by 
the occurrence of pairing, thereby making the order 
that should prevail in the low-temperature under- 
doped region rather elusive. Besides the possibility 
of unusual ordered states [2,5], experiments indicate 
that charge ordering (CO), possibly in the form of 
stripe or checkerboard textures, is indeed present in 
underdoped superconducting cuprates. A survey is 
provided in Refs. [1,6]. 

Based on the experimental evidence and on theo- 
retical results within models for strongly correlated 
electrons, a scenario was proposed [4], where quan- 
tum criticality near optimal doping and CO in the 
underdoped region merge in a unified scenario. Fur- 
thermore, due to the closeness to an antiferromag- 
netic (AF) Mott-insulating phase at very low dop- 
ing, also the presence of substantial AF spin fluc- 



tuations could markedly affect the physical proper- 
ties of these systems [7] . The roles of AF and charge 
fluctuations are not necessarily exclusive: The for- 
mer could rather be emphasized by CO, and be 
"enslaved" in the charge-poor regions produced by 
the modulation of the charge profile, thus extending 
their relevance to doping levels substantially higher 
than the AF onset. 

Within a model of strongly interacting electrons 
coupled to a dispersionless phonon [4], for realistic 
model parameters, the CO instability first occurs 
along the (1,0) or (0,1) directions (i.e., along the 
Cu-0 bonds) with typical wavelength of a few lattice 
spacings. Accordingly, near the QCP the fermion 
quasiparticles (QP) are coupled to low-energy CO 
collective modes (CM), and are strongly affected 
near the "hot" spots at the Fermi surface, which 
are mutually connected by th critical wavevector q c 
and reside near the points (±7r, 0) and (0, ±7r) of the 
Brillouin zone. A similar phenomenology is found 
in the case of AF spin fluctuations, although with a 
different characteristic wavevector, q s « (7r, tt). The 
effects of the exchange of nearly critical CM have 
been extensively studied within phenomenological 
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models with QP-CM coupling. In this way, specific 
features in ARPES (e.g., the QP self-energy [8], the 
kink [9] and the isotopic dependence of the disper- 
sions [10]), and anomalous isotopic effects [11] have 
been coherently explained within the CO-CM sce- 
nario. 

In this paper, we explore the optical and Raman 
signatures of the CO and of the (enslaved) AF spin 
fluctuations, in particular (a) the possibility to de- 
tect spectroscopic features displaying scaling prop- 
erties as fingerprints of quantum critical CM and (b) 
how the symmetry properties of the Raman selec- 
tion rules determine the relevance of the spin and/or 
charge CM at different energy scales. 

2. Critical collective modes in optical spectra 

Near the CO instability (i.e., near the QCP at 
optimal doping and in the underdoped, on the verge 
of CO), the CO CM have a small mass and provide 
a "cheap" reservoir of excitations, which can easily 
affect the various spectroscopic probes. 

Recently, we performed a microscopic calculation 
of the optical conductivity [12], based on a 

conserving approximation of the current-current re- 
sponse, which is crucial to correctly describe opti- 
cal absorption. It was important to realize that for 
a parabolic electron dispersion and in the absence 
of quenched impurities (clean case), the CM cannot 
be purely electronic and must involve, e.g., phonons, 
otherwise the conservation of electron momentum 
would entail a vanishing a(ui) at w > 0. Phonons 
modify the purely electronic relaxational form of the 
CO CM propagator, yielding an effective interaction 

z/(q - c\ c ) z + \u n \ + uj^/il + m 

among the QP. Here, u n are boson Matsubara fre- 
quencies, g is the QP-CM coupling, v is a fermion 
scale, the high-frequency cutoff relates fermion 
and phonon scales [12] and encodes the dynamical 
nature of the phonons, and m cx is the CO CM 
mass, proportional to the inverse square correlation 
length measuring the distance to the CO transition. 
A similar effective interaction describes also the AF 
CM, with a proper redefinition of the parameters, 
and yields very similar results for <t(uj). For Q < oo, 
a finite response was found, with various regimes for 
cr(co) [12]. At low temperature T, the near absence 
of dissipation for ui < ft strongly suppresses a(u)), 
yielding a peak which is fixed at uj ~ f2 and does 
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Fig. 1. Diagrams for the CM corrections to the response 
function. The full dots represent the current vertices in the 
case of optical conductivity or the Raman vertices. The solid 
lines represent the QP propagator (with or without impurity 
scattering) and the dashed lines represent the interaction 
mediated by CO or AF CM. Self-energy, vertex, and Aslama- 
zov-Larkin-like contributions are labeled by S, V, and AL, 
respectively. 

not scale with m when m < f2. On the other hand, 
if m > Q, a peak related to the excitation of CO 
CM appears at uj ~ m. At higher T, the thermal 
excitation of the diffusive CM rapidly fills the low- 
frequency spectrum, embedding the finite- frequency 
peaks into a broader peak at uj = 0. Therefore, the 
scaling behavior may be masked by non-universal 
effects, and we concluded that the lack of scaling 
peaks in optical spectra cannot be taken as evidence 
against quantum criticality. 

Here, we address the question whether and how 
the presence of impurities (dirty case) modifies the 
above conclusions. We determine the conditions for 
a finite-frequency peak to be visible and display scal- 
ing via its mass-dependent position, when an impu- 
rity scattering rate T broadens the QP propagator in 
Fig. 1. In these diagrams the impurity vertex correc- 
tions cancel out, while, owing to the finite q c , no dif- 
fusive impurity ladders can be inserted. We assume 
a CM mass m oc T, as appropriate for the quantum- 
critical regime, and calculate the corrections Sa(u>) 
arising from the current-current diagrams of Fig. 1. 
We first discuss the purely perturbative result, and 
then introduce a suitable resummation scheme to 
avoid the shortcomings of perturbation theory. 

Impurities enlarge the free-electron 5(tu) response 
into a Drude peak of width T. If T is of the order of 
(or even larger than) m and/or f2, this Drude term 
embeds the peaks due to CM scattering and makes 
the critical behavior of the fluctuations difficult to 
detect. Therefore, a first conclusion is that the ob- 
servation of peaks due to CM is only possible when 
disorder is sufficiently weak. To further explore the 
effects of critical CM in optics we accordingly in- 
troduce a rather small disorder scattering rate T < 
m, fL Disorder produces a finite negative part in the 
CM perturbative correction, 5a(ui), at low uj. This 
is naturally so, because the conservation of the op- 
tical weight implies that 5a has zero total weight, 
/ duiSa(uj) = 0. Therefore, when the dirty QP are 
dressed by CO CM, the spectral weight added at 
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frequencies ui > T has to be subtracted at u> < T. 
This feature was already present in the clean case, 
where, however, the negative part was 5(w)-like and 
did not affect the spectra at finite 10 [12]. 

In the presence of impurities, the negative part 
of 8a may prevail despite the Drude peak even for 
small or moderate QP-CM coupling g, yielding an 
unphysical negative absorption at low lu. This arti- 
fact can be cured by an appropriate resummation 
scheme, which we propose here. Specifically, we con- 
sider a memory function [13], which treats on equal 
footing the QP scattering due to impurities and to 
nearly critical CM, M(lo) = iT - w^f(w) /W 
where Xcm™ * s the clean CM correction to the 
current-current response (see Fig. 1) and W is 
the free QP optical weight. The corresponding re- 
summed response is 



Wlj 



M{u>y 



(2) 



which reproduces the perturbative result at high oj, 
where XclT / w < In the absence of QP-CM cou- 
pling Xcm™ = an d we recover the bare Drude re- 
sponse xd — Wu>/ (oj + iT) . Although our resumma- 
tion scheme is no longer conserving, it does not spoil 
the previous conclusions on the conditions for opti- 
cal peaks scaling with to. Figs. 2(a) and (b) report 
the results for the resummed a(uo) = -Imx 3 j(w)/w 
for the cases to > f2 and m < Q, respectively. 

Our analysis shows that, in the presence of a small 
amount of impurities, the spectra broaden, but the 
generic feature persists, and the absence of a In- 
dependence in the peak position [see Fig. 2(b)] can- 
not be used to infer the absence of critical CM (with 
quantum-critical mass m oc T). The CM become 
fully visible in optical spectra only at low T and 
when the CM mass is larger than a typical non- 
universal dissipation scale Vt, [see Fig. 2(a)]. Note 
that the displacement of weight from oj < T to 
higher frequencies, produced by 8a , can result in a 
narrowing of the Drude peak. In the analysis of ex- 
periments, this could lead to an erroneous overesti- 
mate of the QP lifetime. The determination of the 
precise amount of narrowing and of the extent to 
which the negative part of 8a could "carve" a dip in 
the spectrum, thereby producing a finite-frequency 
peak despite the rapid thermal filling mentioned 
above, requires a reliable conserving resummation 
scheme to treat the low-frequency absorption, which 
is presently not available, due to the lack of a small 
expansion parameter. 
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Fig. 2. Optical conductivity for a system of QP coupled to 
CO CM modes in the presence of finite impurity scattering 
rate T = 3 cm" 1 , at various temperatures and CM masses, 
for (a) Ti = 30 cm" 1 , and (b) Ti = 330 cm" 1 . 

In summary, within our approach a sufficiently 
strong disorder may effectively wash out the effects 
of CM, thereby hiding criticality. When disorder is 
weak enough to leave the peaks at low T visible, im- 
purities do not qualitatively modify the findings of 
the clean case: The key role in the non-universal be- 
havior of a(uj) is played by charge conservation in 
the presence of an extrinsic (e.g., of phonon origin) 
energy scale for dissipation fl. The value of the crit- 
ical scale to in comparison to il then determines the 
presence or the absence of a critical T dependence 
of the optical features. In cuprates we estimate O ~ 
100 cm" 1 [12]. 

3. Critical collective modes in Raman 
spectra 

The natural question arises on the role of dissi- 
pation on the behavior of other spectral quantities, 
which are not "protected" by charge conservation. 
Within this respect, Raman spectroscopy provides 
a nice playground with the additional advantage of 
a momentum-dependent filtering of the spectra. A 
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key difference with respect to the optical conduc- 
tivity is that, properly adjusting the polarization 
of the incoming and outgoing photons, the Raman 
response involves form factors selecting specific re- 
gions of the QP Brillouin zone [16]. This helps to 
distinguish nearly critical AF and CO fluctuations 
in Raman spectra, which would look very similar in 
optical spectra [14]. 

To compare Raman and optical spectra, we con- 
sider here the contribution to the Raman response 
due to AF and CO CM. Since we have shown that 
impurity scattering should be weak enough not to 
hide the effects of CM in optics, we henceforth focus 
on the role of CM and neglect impurities. Specifi- 
cally, we consider the clean Raman response Im% in 
the Big and Big channels, where photons couple to 
QP via the vertices 7^ = cos k x —cos k y and jb 2 = 
sin k x sin k y . While jB lg is large (and therefore se- 
lects excitations) in the antinodal regions around 
the (±7r, 0) and (0, ±7r) points of the Brillouin zone 
and changes sign in the nodal regions k x = ±k y , 
r )B 2g is large around the nodes and changes sign in 
the antinodal regions. These symmetry properties 
induce leading-order cancellations between the Ra- 
man response diagrams depending on the momen- 
tum exchanged by the CM. For instance, we find 
that the S and V diagrams (see Fig. 1) cancel each 
other in the B\ g channel, when typical momenta q c 
of CO excitations are exchanged, because 7^ has 
the same value at the two hot spots connected by q c , 
and the S and V diagrams have intrinsically oppo- 
site sig rQ. In this case, the only relevant contribu- 
tion arises from the AL diagrams, which were con- 
sidered in Ref. [15], to explain the Raman absorp- 
tion anomalies observed in La2- ;c Sr :r Cu04 (LSCO) 
with x = 0.02,0.10 [18]. On the other hand, AF 
CM with typical wavevector q s ss (n, tt) give qual- 
itatively wrong results. In Ref. [15], a diffusive CM 
propagator was considered, mediating an effective 
interaction similar to that in Eq. (1), but with a 
different high-frequency behavior. Here, having in 
mind a comparison with the results on cf(lu), ob- 
tained using Eq. (1), we reproduce the analysis of 
Ref. [15] using the high-frequency cutoff fi. As it 
can be seen in Fig. 3(a), the experimental data are 
well reproduced setting fl = 200 cm -1 and using m 
as the only fitting parameter, once the overall in- 
tensity is set by the curve at the highest tempera- 



1 This cancellation of S and V diagrams is reminiscent of the 
loop cancellation occurring in the density-density response 
of systems with strong forward scattering [17]. 
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Fig. 3. (a) Raman absorption in the B\ g channel for a 
clean system of QP coupled to CO CM at temperatures 
(from top to bottom) T=36,65,94,124,188 K. The mass m(T) 
needed to fit the data is reported in the inset. A dissipa- 
tion scale f! = 200 cm -1 is used, (b) Same as in (a) with 
T=36,94,124,188 K, but with a smaller dissipation scale, 
Ti = 20 cm- 1 . 

ture. Again, m(T) behaves as expected for a system 
crossing over from the quantum critical region to 
the low-T nearly ordered phase. Clearly, due to the 
inherently non-conserving character of Raman ex- 
citations, the scaling behavior of the peak is visible 
also when m < f2. For the sake of completeness, we 
also discuss the case m > f2, although it does not 
apply to LSCO. In Fig. 3(b) we report the B\ g Ra- 
man spectra (again coming from the AL diagrams 
only, for symmetry reasons) for the same masses as 
in Fig. 3(a), but a smaller fl = 20 cm -1 . Here, a 
narrow peak appears at a frequency uj ~ Q, above 
which the CM changes from diffusive to propagating 
[see Eq. (1)]. The mass energy scale m is still visible 
as a shoulder (for low and intermediate T) or as a 
broad peak (at the largest temperature T ~ 188 K). 

4. Symmetry properties of Raman spectra 

As mentioned in the previous section, the specific 
symmetry of Raman vertices was exploited in Ref. 
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AL^O, S + V&O 


AL ss 0, S + 


AF CM 


AL&0, S + V^O 


AL fa 0, 5 + V ss 



Table 1 

Leading corrections to the Raman absorption due to the 
processes represented by the diagrams of Fig. 1, in the B\ g 
and B2 g symmetries. 

[15] to explain the anomalies observed in LSCO [18] . 
In this section, we focus on those features of Ra- 
man spectra which can be discussed based only on 
symmetry considerations. These provide rather ro- 
bust arguments on the relevance or irrelevance of the 
leading contributions from CO or AF CM in the var- 
ious channels when the CM are nearly critical and 
strongly peaked at their characteristic wavevectors. 
On the other hand, far away from criticality, when 
m is large and the CM mediate smooth and weak in- 
teractions, the symmetry arguments are of little use 
because the CM correction to the response is small 
and substantial contributions come from regions of 
the Fermi surface away from the hot spots. 

By direct inspection of the symmetry properties 
of the QP loops entering the AL diagrams (see Fig. 
1), one realizes that the leading contributions of AF 
CM with q s = (n, ir) are non-zero only in the -fA lg — 
cos(k x ) + cos(k y ) symmetry [19]. In the B ig channel, 
a leading contribution to Raman absorption from 
the AL diagrams can only arise if the CO CM have 
a characteristic wavevector q c in the (1, 0) or (0, 1) 
directions [15]. As far as the S and V diagrams are 
concerned, as anticipated in the previous section, 
the leading contributions in the B ig symmetry can- 
cel when the CO CM is considered. In the Bi g sym- 
metry, instead, the S and V diagrams sum up (since 
now the Raman vertex has opposite sign in the two 
hot spots connected by q c , avoiding loop cancella- 
tion) and provide a finite leading contribution. The 
opposite situation is found when the CM wavevector 
is q s « (7r, 7r), typical of AF CM near the AF insta- 
bility: due to the symmetry of the form factors, the 
S and V diagrams sum up in the B lg channel, while 
they cancel at leading order in the B^ g channel. The 
situation is summarized in Table 1. 

Since the leading S+V contributions diverge as 
u) — > 0, a resummation is required to obtain a fi- 
nite Raman response. We adopt the same memory- 
function scheme of Eq. (2) (here we consider the 
clean case T = 0). We calculate the corrections to 
Raman absorption due to both AF and CO CM in 
the two channels. Specifically, in the inset of Fig. 
4(a) we report the changes in the B lg and B 2g Ra- 



man absorption due to a progressive reduction of 
the CM mass, keeping T and QP-CM coupling fixed. 
This situation is representative of experiments com- 
paring overdoped samples (large CM mass) with 
optimally and underdoped samples (with progres- 
sively smaller CM mass). The main effect is a shift 
of weight from low to high uj. The same shift oc- 
curs upon reducing both T and CM masses, as can 
be seen in the main panel of Fig. 4(a), representa- 
tive of the quantum critical regime for the CM. No- 
tice that, according to the standard expression of 
Raman vertices [16], in this figure the intensity of 
the Bi g spectra has been rescaled by a factor t' /t ~ 
0.2, where t' and t are the next nearest neighbor 
and the nearest neighbor hopping parameters of a 
tight-binding model, respectively. Moreover, while 
the £?2 3 channel never acquires leading AL correc- 
tions from AF or CO diagrams (see Table 1), the 
Bi g absorption may be enhanced if CO CM become 
increasingly critical when the doping is reduced. In 
this case, their additional AL contribution appears 
at frequencies u> < £1 ~ 200. . .300 cm -1 . In Fig. 
4(b) we report the B\ g absorption arising from the 
superposition of AF CM acting via S+V processes 
(resummcd within a memory-function scheme to re- 
move their singular behavior at low frequency) and 
CO CM in the AL channel. Notice that here CO and 
AF CM have the same mass to schematize their en- 
slaved quantum critical behavior. 

Several of the above features are indeed gener- 
ically observed in LSCO and Bi2SrCa2Cu20s 
(BSCCO) [20]. In the B\ g channel, upon going from 
overdoped to underdoped systems, weight is shifted 
to higher frequencies. According to our scheme, 
this effect could be attributed to AF CM acting via 
S+V processes, and becomes progressively more 
important upon undcrdoping, both in LSCO and 
BSCCO. The effect of CO CM via S+V processes is 
nearly cancelled in the B\ g channel and the anoma- 
lous absorption, which we interpret as due to the AL 
contribution, is rather weak even around optimal 
doping. Upon entering the underdoped region, the 
marked tendency of LSCO to form stripes makes 
the AL contribution of CO to the B\ g channel 
more important at low frequencies (uj < 200 cm -1 ). 
Then, while the B\ g absorption stays moderately 
suppressed by the AF S+V contributions at inter- 
mediate frequencies, it acquires a nearly critical 
additional contribution from the AL processes [15] 
and displays the anomalous low-frequency absorp- 
tion shown in Fig. 4(b). Notice that a similar addi- 
tional anomalous absorption is absent in BSCCO 
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Fig. 4. (a) Raman absorption for a clean system of QP 
due to S+V processes. Solid lines: B\ g channel with QP 
coupled to AF CM with Oaf = 500 cm -1 ; dashed line: Big 
channel with QP coupled to CO CM with Q = 200 cm -1 . 
AF and CO CM have masses ms and mc, respectively. Inset: 
Same as in the main panel, but with fixed T = 120 cm -1 
and two different masses of the CM: m = 70 cm -1 (thick 
lines) and m = 200 cm -1 (thin lines), (b) Total B\ g Raman 
absorpion due to AF CM (S+V diagrams) and to CO CM 
(AL diagrams, also shown separately as dashed lines). Masses 
and temperatures are the same as in panel (a). 

samples. We argue that in BSCCO samples the ten- 
dency to form nearly static stripes is much weaker, 
so that the low-frequency AL absorption could be 
masked by other processes competing with CO (e.g., 
pairing without superconducting phase coherence). 
Another possibility (not necessarily exclusive) is 
that in BSCCO the CO texture has the different 
form of a checkerboard (previously named "egg- 
box" in the literature [21]). In this case, symmetry 
arguments can be found showing that checkerboard 
charge fluctuations should couple much less to the 
QP loop entering the AL diagrams of Fig. 1 even for 
the B\g symmetry. Of course these non-trivial sig- 
natures naturally require more specific parameter 
tuning and more realistic model calculations (e.g., 
with specific band structures) and treatment of the 
cutoffs to improve the quantitative agreement with 
the behavior of various samples. 

In conclusion, the robust symmetry arguments 



outlined above account for the generic behavior of 
the Raman spectra at different dopings and in dif- 
ferent channels and, contrary to the optical conduc- 
tivity, allow to distinguish the effects of charge and 
spin modes. 
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